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Abstract 

We show how to make 0(a) corrections in the bare parameters of 3-D lattice SU(2) xU(l) 
Higgs theory which remove 0(a) errors in the match between the infrared behavior and 
the infrared behavior of the continuum theory. The corrections substantially improve the 
convergence of lattice data to a small a limit. 
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1 Introduction 



It is believed that the baryon number of the universe may have originated at the electroweak 
phase transition. The details of the mechanism depend substantially on the details of the 
phase transition itself, so for several years there has been growing interest in a quantitative 
understanding of the electroweak phase transition. Ordinary perturbation theory |l] has 
proven useful only when the Higgs mass m# <C mw the W mass, and outside of this regime 
the only reliable methods are nonperturbative, such as lattice Monte-Carlo techniques. 

The reason for the failure of perturbation theory was elucidated by the work of Farakos 
et. al. 0, |3|, M, who have shown how, to very good approximation, the thermodynam- 
ics of the standard model near the electroweak phase transition can be described using a 
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three dimensional path integral for purely bosonic SU(2)xU(l) Higgs theory. The theory is 
super-renormalizable, and hence its ultraviolet is well behaved, but its infrared is potentially 
nonperturbative, leading to the breakdown of perturbation theory as a tool for exploring 
the basically infrared physics of the phase transition. However, this same feature makes 
the numerical investigation of the nonperturbative physics particularly feasible; fine lattices 
are possible because the system is low dimensional, the couplings fall quickly with lattice 
spacing, and there is no need for the (numerically expensive) inclusion of fermionic degrees 
of freedom. 

Several groups have conducted numerical experiments on the three dimensional system 
to study the phase transition || || |7|, || [lOfl . Each study used the Wilson action and 
the bare relations between lattice couplings and wave functions and the 3-D continuum 
parameters. In each case there were substantial linear in a (lattice spacing) errors which 
had to be removed by extrapolation over several values of a. The extrapolation process is 
numerically expensive because it demands very fine lattices, and the errors in the final result 
are dominated by 0(a 2 ) errors from the coarsest lattice used, and statistical errors in the 
finest lattices. It would be better to have an analytic understanding of the origin of the 
linear in a corrections, which would make it possible to prevent them. One could then take 
data only on lattices of intermediate coarseness, obtaining smaller 0(a 2 ) errors and better 
statistics with less numerical effort. This paper presents an analysis of these 0(a) errors. 

The 0(a) errors in measured infrared phenomena arise because the infrared theory is 
clothed by interactions with more ultraviolet modes, but the clothing differs between lattice 
and continuum theories, both because the most ultraviolet degrees of freedom are absent on 
the lattice, and because the ultraviolet lattice modes have incorrect dispersion relations and 
extra (lattice artefact) interactions. Since the dispersion relations and couplings of the lat- 
tice theory agree with those of the continuum theory in the infrared, the difference between 
clothing in the two theories is all ultraviolet, and therefore perturbatively computable. Also, 
since the origin of the difference is ultraviolet, it should be expressible as an operator ex- 
pansion, and at the desired level of accuracy only the superrenormalizable terms are needed; 
these can be compensated for by a shift in the bare parameters of the theory. 

Dimension two operators, such as the mass squared and the $ 2 operator insertion, receive 
linearly divergent corrections at one loop; logarithmic corrections, computed in O, at 
two loops; and 0(a) corrections at three loops. The problem of determining these 0(a) 
corrections is difficult. However, if one is only interested in properties of the phase transition 
such as jumps in order parameters, the latent heat, the surface tension, etc. then it is not 
necessary to know these to high accuracy. To compute the jump in $ 2 , for instance, one 
only needs to account for the multiplicative correction to the operator insertion, not the 
additive correction, which is common to the two phases. Similarly, precise knowledge of 
m 2 H is only necessary to compute very precisely the equilibrium temperature. We will not 
pursue improvement in the determination of these quantities here, but will concentrate on 
improving the precision of the measurement of quantities related to the strength of the phase 
transition, which is more physically interesting. 

The outline of the paper is as follows. In Section |2| we present the problem, establish 
notation, and give the results. Readers uninterested in the details of the calculation can then 
skip Section |3], which enumerates the required Feynman diagrams, presents the calculation 
of a few, and gives results for the others. The conclusion, Section f|, contains an example 
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of using the improvement technique to reinterpret existing data. The first two appendicies 
contain integrals, identities, and Feynman rules needed in Section [3], and the last presents 
the results for the abelian Higgs model, and for the standard model but including the adjoint 
scalar (A ) fields. 

2 Set-up and Results 

The dimensional reduction program shows that the infrared thermodynamics of the standard 
model are described up to small error by the partition function 

Z = J V^VAiVB ie xp(-l3 J d 3 xH), (1) 
H = ^(^ + —^7^ (2) 

Here A4 and £>« are the gauge fields and and are the field strengths of the SU(2) 
and U(l) gauge groups respectively, and the couplings and wave function normalizations are 
determined with respect to physical measurables as discussed in 0. Only m 2 H depends on 
the (MS) renormalization point fi. 

The corresponding (Wilson) lattice action, in lattice units, is 



Z 



J VUVuV® ex P (-(3 L ]T H L ) , (3) 



Hl = ^E(l-^TrP tJ ) + -J^5:(l-Re^) 

i<j 1 tan u w i<j 
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+Z„ PT($a(^) - <S> a {x)U iab {x)u ibc {x)<S> c {x + z)) + -^§1 



(Ax - 5X L )Zl 
4 



{KY , (4) 



where $ is treated as four independent real entries $ a , U and u are the SU(2) and U(l) 
link variables respectively, and and are the SU(2) and U(l) lxl plaquettes. The 
notation for <3> a is equivalent to the common notation in which the Higgs field is written as 
a magnitude times an SU(2) matrix, and each sum on a is replaced by (l/2)Tr. 

The (naive) tree relations between the lattice parameters and the 3-D continuum ones 

are 

H = If, (5) 
9*0, 

Z A = 1 (6) 
Z B = 1 (7) 
Z® = 1 (8) 

Al = ^, 5\ L = 0, (9) 

g 2 
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42 = 8^£ = 4^, 



where is the $ 2 operator insertion. 

The super-renormalizability of the theory ensures that, in the limit a —> (or (3^ — > 
oo), using these relations-except for Eqs. (|T0| ) and QTT|), which must be corrected for one 
and two loop divergences-will make the thermodynamics of the lattice system match the 
thermodynamics of the continuum theory. However, the approach to the thermodynamics 
of the continuum theory need not be very fast; while we are assured that it is power-law, it 
proves to be 0{Pj}) (or 0(a)). This is because one loop effects introduce 0((3~[ l ) corrections 
to the righthand sides of Eqs. (H) through ([□]), as discussed below. In the case of Eqs. ( |TOl) 
and (|Tl|) , the 0((3j}) correction swamps the 0(P£ 2 ) tree value; there is a linear divergence, 
which should be tamed in order to control the lattice theory. This correction, and 2 loop 
0(/?^ 2 ln/? £ ) corrections, are computed in [§, [□]]. The corrections to the other equations 
have not previously been computed. 

To see how O^^ 1 ) corrections arise in the infrared behavior, consider some Feynman 
diagram needed in the evaluation of some infrared quantity of interest, such as the effective 
potential or the domain wall surface tension. The dominant contribution arises when all 
loop momenta are on order the nonperturbative scale p ~ g 2 T, or in lattice units p ~ (5j}- 
At these momenta the theories match except for O^J 2 ) nonrenormalizable operators, which 
are not important at the desired level of accuracy. There are also contributions where one or 
more loop momenta are "hard," p ~ 1 in lattice units. At such momenta the lattice and con- 
tinuum theories differ significantly. However, at these momenta the theory is perturbative, 
by the same power counting argument which establishes its super-renormalizability. For the 
wave function normalizations of Eq. (|4]), each vertex contributes a power of (3l and each 
propagator contributes a power of so the result is suppressed parametrically by one 
power of 1 for each hard loop. The value of a hard loop with all soft incoming momenta 
can be computed once and for all in an expansion in the momenta of the incoming soft lines, 
so the difference between the lattice and continuum hard loops can be written as an oper- 
ator product expansion. Only dimension 2 and 4 (super-renormalizable) operators will be 
important at 0(f3~[ l ) for the evaluation of the remaining soft diagram, as higher dimension 
operators just contribute small corrections to the nonrenomalizable operators which could 
already be neglected. And the super-renormalizable operators can be absorbed by O^j}) 
counterterms in the tree level Lagrangian, making the remaining soft integrations match 
between the lattice and continuum theories, up to O^fiJ 2 ) corrections. In other words, we 
resum the class of all diagrams with non-overlapping hard loops by collapsing the hard loops 
into operator insertions and applying counterterms in the Lagrangian to absorb the differ- 
ence between the lattice and continuum values. The lattice and continuum perturbation 
series then match up to 0(/3^ 2 ) corrections when we sum over the subset of diagrams with 
isolated hard loops, rather than just those diagrams without hard loops, which is the subset 
for which tree relations between lattice and continuum parameters lead to matching contri- 
butions. Since overlapping hard loops are suppressed by two powers of by the same 
power counting used earlier, the infrared theories will now match at 0(/3 i ^ 1 ). 
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An alternate way of phrasing the same argument is to consider how the lattice theory 
can arise from integrating out degrees of freedom from the continuum theory. We want the 
infrared behavior to match when we integrate out these degrees of freedom. At tree level 
one just relates the couplings and wave functions by naive dimensional analysis. But to 
achieve higher accuracy one must take into account the loop effects which were provided by 
the degrees of freedom no longer present, and the change in loop effects due to the inevitable 
changes in the theory at the lattice scale such as the new lattice dispersion relations and 
added (lattice artefact) couplings. The influence of these loop effects on the infrared physics 
should be computed; as discussed above it can be expressed at the desired level of accuracy 
as perturbatively small operator insertions, of which only the super-renormalizable ones are 
relevant. One should then shift the values of the lattice theory Lagrangian parameters to 
absorb the difference between the lattice and continuum contributions, ie one replaces loop 
effects which would appear in the continuum theory but which are absent in the lattice 
theory with shifts in the lattice parameters which reproduce these effects. The corrections 
are perturbative with perturbation expansion parameter fi^ 1 , because they all arise from 
momentum scales on order or greater than the reciprocal lattice spacing, and the power 
counting discussed above applies. (The only exception would be if there were infrared diver- 
gences, but since the infrared lattice and continuum theories match, the difference in loop 
contributions between the two theories will be free of infrared divergences.) At the desired 
level of accuracy we need to perform a one loop match. 

A complication is dimension 2 operators. The most general renormalization of these 
operators is 



'/fX,renormalizcd 



op,renormalized 

In addition to a multiplicative renormalization, Z m and Z op , there is an additive renormal- 
ization, 5m 2 HL and 5&l p . The problem is that in lattice units the tree level m 2 HL , and the 
expected value of <3> 2 due to the interesting infrared physics, is 0(P^ 2 ). The O^j}) correc- 
tions 5m 2 HL and 5$ 2 p will swamp these. That is why it was considered necessary to carry the 
calculation of these counterterms to two loops. To determine the counterterms to O^fi^ 1 ) 
accuracy relative to the tree value would require a prohibitively hard three loop calculation. 
However, in practice we only need to know 5m 2 HL with precision if we want to determine 
the phase transition temperature with high accuracy, since the relation between the 3 di- 
mensional continuum value of m 2 H and the physical, 4 dimensional Higgs mass squared is 
strongly temperature dependent. In practice one tunes m 2 HL ren — 5m 2 HL to find the phase 
transition temperature; to find the absolute temperature this represents, one must know 
5m hl with precision, but to determine the splitting between the equilibrium temperature 
and the spinodal temperatures, for instance, one only needs to know Z m ; only this quantity 
is relevant for details of the strength of the phase transition. Similarly, we need to know 
5Q 2 to find the absolute value of 2 in one phase; but the question which is more interesting 
physically is the difference in 2 between the two electroweak phases, and 5$ 2 p cancels out 
in such differences, and only Z ap is relevant. Furthermore, because the mass squared and 



lQ(3 2 m 2 H 
4/3 2 2 



+ Sm 2 HL , 



op 



J op 



9 2 Pl 



5$ 



op 



(12) 
(13) 
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Figure 1: Relation between bare and clothed propagator at one loop. 

the <J> 2 operator insertion couple to the same operator, Z m = Z op ; so henceforth we will 
only consider Z op (though we will also reproduce the (known) one loop value for 8m 2 HL as a 
natural biproduct of our renormalization calculation). 

Now we proceed to set up the computation of the renormalization constants needed in the 
match between continuum and lattice theories. We begin with the scalar propagator. In the 
continuum theory the value of the scalar propagator is G(p 2 ) = l/(p 2 + m 2 H — 7r^ t c{p 2 )), with 
7r$ c* the continuum theory self-energy computed using the clothed vertices and propagators, 
expressed in lattice units. In the lattice theory the bare propagator is modified by Z J , so 
the clothed propagator is (see Figure [I]) 

G = 2 ■ p- . 2 (1 + ^ L G) , (14) 

V + m HL ~ °™>HL 

and hence 

(Z& - l)p 2 + (Zj, - l)m 2 HL - Z$5m 2 HL = 7r #)i - 7r$ iC . (15) 

Hence one can determine 8m 2 HL and Z$ from the difference between the lattice and continuum 
self-energies at O(p ) and 0{p 2 ). The self-energies should be computed using clothed vertices 
and wave functions, which coincide (by hypothesis) in the infrared; so the difference will be 
ultraviolet dominated, and here the clothed vertices and wave functions can be replaced by 
bare ones (without Z and 5\ corrections) with an 0{j3j j l ) error, which would be accounted 
for in a full two loop calculation and will only lead to an 0(/5^ 2 ) error in and can 
therefore be neglected. Similarly, the (Zq> — l)m 2 HL can be dropped {m 2 HL being 0(/?^ 2 )), 
and the Z$> in front of Sm 2 HL can be set to 1. 

Next, denoting the loop contribution to the amputated 1PI four point function at small 
(zero) external momentum as —V\ (the minus sign because the contribution from the scalar 
self-coupling is —A), and again denoting the lattice and continuum values with L and C 
subscripts respectively, one demands that the actual strength of the scalar self coupling 
coincide between lattice and continuum theories, 

(X L - 8\ L )Z% + V X , L = X L + V x ,c , (16) 

and hence, at leading order in /3~£ , 

SXl = V\,l - V x , c + 2X(Z^ - 1) . (17) 
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Z A P 2 + 71 A 



Full effect goes as (Z + V A Y 



1 

Z A v 2 +n A 



z^+v A 

Figure 2: The full influence of a gauge particle propagating between two lines depends on 
the vertex and propagator renormalizations. Here and throughout SU(2) lines will be curly 
and U(l) lines will be wiggly. 

Again, the V\ are computed with the clothed vertices and propagators, but can be computed 
using the naive bare propagators and vertices because the difference will be ultraviolet dom- 
inated. 

To find Za and Zb one must compute not only the gauge field self-energy corrections -ka 
and 7Tb, but also the corrections to the gauge-scalar vertex, which we will denote as Va, Vb- 
(It does not matter which vertex is considered, by gauge invariance; we choose the scalar 
vertex because it is the easiest to compute.) Za should be chosen so that the complete effect 
of a gauge particle propagating between two lines, say two scalar lines, is the same between 
the lattice and continuum theories, see Figure 0. The strength of the vertex at each end of 
the propagator is multiplied by Z$ + Va,l in the lattice theory and 1 + Va,c m the continuum 
theory; the Z$ factor arises because its appearance in the action modifies both the wave 
function and the coupling to the gauge field. The modification of the propagator is similar 
to that for the scalar, and we find that to match the full modification (a vertex at each end 
and a propagator in the middle) we must choose Za such that 

(gg + Va,l? 1 + Va.c , 1a , 

~V 2 ~ ~2 ' 

or, at leading order in fij} , 

Z A -1 = 2{Z* - 1) + 2{V a ,l - V A ,c) + • (19) 

The equation for Zb is analogous, and in this case the vertex and Higgs wave function 
contributions will turn out to cancel. 
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The $ 2 operator multiplicatively renormalizes due to the diagrams shown in Figure || 
because in every diagram where the operator appears on an infrared Higgs field line, there 
is a corresponding diagram where the $ 2 insertion on that line is replaced with the one loop 
insertion shown in that figure. If the sum of the contributions of the diagrams there, at small 
(zero) external momentum, are denoted by L$,z,, and the continuum contribution is L$,c, 
then the renormalization of the $ 2 operator is 

Zop — 1 — {L<s>,l — L$ t c) j (20) 

where each L should be computed with clothed propagators, vertices, and operator insertions. 
Again, at the correct value of Zqp the infrared definition of the operator insertions will 
coincide, so there is no infrared divergent contribution to the difference; the difference is UV 
dominated and can be computed, at the desired level of accuracy, with the bare vertices and 
propagators. 

It is possible to scale $ by a constant k without changing the physics, provided one 
multiplies Z$ by k~ 2 and Zqp by /t -2 . Only the combination Z$,Zq\> must be meaningful 
and gauge invariant. It is convenient to perform such a scaling of <3> to set Zqp to 1, in which 
case one should add — £$,c) to Z<$, at the very end of the calculation. 

The diagrams needed for the renormalization are presented in the next section, and their 
values are tabulated there. All but two are performed in lattice Lorentz gauge with a general 
gauge parameter, to check gauge parameter independence; the exceptions are the gauge field 
self-energy diagrams with all gauge vertices, which are only performed in Feynman gauge; 
so we have not yet checked the gauge invariance of Za- The results of the calculation are 

z A _ 1 = K i (11± + i = + IV (21 ) 

HL \ 3 4vr 12 4tt 3 J K J 

Z B -1 = ^ tan 2 9. (-1^ + 1^ + 1), (22) 

Z^Zol-l = PE 1 ((9 + 3 tan 2 6^ - 6A) ±- + 3 + ta f ® W , (23) 
= /5Z 1 (6 + 2tan 2 6 w + 6A L )^, (24) 

47T 

S\ L = pi 1 ^(-12A| -4-2(1 +tan 2 9vi/) 2 + 6A L (3 + tan 2 evi/)) ^+ (25) 

3 + tan 2 6 Wx S\ 
3 XL ^J ■ 

The constants S = 3.175911536 and ^ = 0.152859325 have the same meaning as in 0, and 
arise from integrals presented in Appendix A. If one is interested in the theory without the 
U(l) gauge group, set tan 2 Q w = and disregard Z B . Equations (pl|)-(|25|) are the main 
result of this paper. 

By applying these corrections to the Hamiltonian of a lattice simulation one should be 
able to remove all 0(a) errors (except in the absolute determination of the order parameter 
and the phase transition temperature, as discussed above). If one already has numerical 
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data, then the 0(a) errors can be removed by re- interpreting the meaning of the answers. 
In the case of SU(2) Higgs theory alone, equating 




naive ' a \ naive) 



(26) 



with 




m HL 0m HL 2 rv2^L,im V V*L (if o ^\ 

z * 2 im p * 4 y^imp) J 



(27) 



one finds as matching conditions 



Pimp Z OP®l 



imp 



imp 



imp 




naive 



(28) 
(29) 
(30) 



One should re-interpret at what temperature the simulation was conducted, what input value 
of Al was used, and what the result for /3|$ 2 wasQ. Note both that the value of \h,imp is 
smaller than the value of \L,naive and that the rescaling of f3 L & 2 reduces its value. Both 
lead to an overestimate (before correction) of the strength of the phase transition, the first 
because the phase transition is stronger at smaller A and the second because the real jump in 
the order parameter is smaller than that using naive conversions. The error in A is the most 
problematic, because at small values of A the strength of the phase transition is strongly A 
dependent (one loop perturbation theory, which applies here, says the jump in $ 2 should 
go as A" 2 ), and at larger values of A the properties of the phase transition may change 
qualitatively as one changes A. Note also that Za — 1 — CI.65/3^ 1 is quite large, and so are 
corrections where Za appears. 

One can get a rough understanding of the magnitudes of Za-, Z B) Z$, and 5m 2 HL by con- 
sidering the tadpole correction scheme of Lepage and Mackenzie jTJJ . Based on a mean field 
theory argument, they propose that the dominant loop contributions at every order can be 
absorbed by making an "educated guess" for the constants Za, etc, as follows: compute the 
mean value of the trace of an elementary plaquette, and call it U (or u for the U(l) plaque- 
tte). Then guess that every term in the action which contains a link should be multiplied by 
J/ -1 / 4 (with the value of U found after the correction has been applied, so U must be found 
self-consistently). Hence, Za should roughly equal U" 1 , Zq> should roughly equal [/ _1 / 4 u~ 1//4 , 
Zb should roughly equal and 5m 2 HL should be about 6(f/~ 1//4 u -1 / 4 — 1) (because the 
hopping term in the Higgs wave function gets corrected, but not the local term). It is quite 
easy to compute at one loop that U — 1 — fij} and u — 1 — tan 2 O w /?^ 1 /3. (This is just 

2 If we include the U(l) subgroup, this procedure does not work quite as well; the value of tan 2 ®\y will 
differ before and after improvement, which might not be desirable. 
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the approximation that the average energy in the plaquette term in the Hamiltonian obeys 
equipartition.) Hence, one guesses that, at lowest order in fi L , 



P L (Z A -1) ~ 1, 



(31) 



p L {Z*Zol-l) ^ - 





(32) 



(33) 




(34) 



In fact, at = 0, the numerical values of the corrections are 



Pl(Z a -1) 

I3l{Z b - 1) 
i5 L {Z^Z \ - 1) 
f3 L 8m 2 HL 



0.65, 

0.350 tan 2 G w , 
0.0786(3 + tan 2 G w ) 
0.505(3 + tan 2 Q w ) , 



(35) 
(36) 
(37) 
(38) 



which are all close. The tadpole argument is least accurate for the SU(2) wave function 
correction Z A — 1. 

The tadpole improvement value for the correction to the scalar self-coupling is that it 
should only change due to the wave function correction. This misses the corrections at even 
powers of Al (though these are all suppressed by the rather small number £/(47r)). These 
corrections are important when Al is fairly small, because results are very Al sensitive in 
this regime. This illustrates a limitation of the tadpole improvement scheme in a theory 
with non-gauge couplings. 

3 Details of the calculation 

Here we will present the calculations of the vertices and self-energies needed in the last 
section in more detail. To illustrate what is involved in the calculation we will present the 
full details of the contribution of loops involving gauge particles to the scalar self-energy and 
of scalar loops to the gauge field self-energy. For all other diagrams we will simply present 
results. 

The diagrams contributing to the Higgs self-energy are presented in Figure |3|. Diagram 
(a) is computed already in |J. It is straightforward and leads to an O(p ) contribution to 
the difference in self-energies of 



Diagram (b) vanishes in the MS regulated continuum theory, so its contribution to the 
difference of self-energies in Landau gauge is 




from (a) . 



(39) 




(40) 
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(a) (b) (c) 




(d) (e) 

Figure 3: Diagrams contributing to the Higgs boson self-energy. The curly line refers to 
the SU(2) gauge particle, the wiggly line to the U(l) gauge particle, and the solid line to the 
scalar. 



where here and throughout p is the momentum of the line, ki = 2sin(fcj/2), and k 2 = 
Z)i=i,2.3 The factor of 3 is a group factor. 

To make expressions more concise, the factor d 3 k/(27r) 3 will henceforth be implied in 
every integral, continuum or lattice. 

At small p, which is the interesting regime, cos(pj) ~ 1 — pj/2, and the contribution from 
diagram (b) becomes 

PL A-*,*] 3 \k 2 Y 2 (k 2 ) 2 J 

The integral of kf /{k ) 2 is just 1/3 the integral of 1/k 2 by cubic invariance, so the contribution 
of the diagram is 

ir*,L-ir* t c = PZ 1 (-6+p 2 )^- from (6). (42) 

The contribution of (d) differs only in counting factors; the 3 is replaced with tan 2 Q w , which 
arises from the U(l) propagator. 

In Landau gauge, diagram (c) on the lattice contributes 

tt* iL from {c)= T UT—lk " • ( 43 ) 

(J L J[-7r,7r] 3 k 2 (p - k) 2 

The continuum contribution is 

n»,c from (c) = - 1 >- . (44) 

Each numerator vanishes at p — 0, so to extract the 0(p 2 ) behavior we need only expand the 
numerator to 0(p 2 ) and use the denominator at p — 0. (Each contribution is then separately 
infrared divergent, but the infrared divergences match and do not matter to the difference.) 
Since 

(2p -k)i = 2 sm( Pi - ki/2) = —ki cos( Pi ) + 2p t cos(/c i /2) , (45) 
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and ki annihilates against the gauge propagator (only because we are in Landau gauge), the 
contributions simplify at small p to 



3 ^ / f cos(ki/2) cos(k j /2)(5 ij - kikj/k 2 ^ 
- from (c) = —2^ 4^ ' ' 



/3l v vi--,-] 3 (Jfe 



2\2 



SR 3 /c 4 



(46) 



Each integral vanishes when i ^ j. When i — j, one may use cos 2 (/cj/2) = 1 — fc 2 /4 and 
cubic invariance to reorganize the contribution as 



3p 2 

tt$,l - 7r*,c from (c) = — 





1 r lr k\ 



3 Vh-,-] 3 (P) 2 ^ 3 fc 4 / 3 J[-*^f k 2 J{-Tr,nf (k 2 f 

"(47) 

The first two integral expressions here are the constants £/47r and S/47T. The last is 
related to these constants by an identity, Eq. (|110|) . The result is 

vr$, L - 7r$,c from (c) = ^- [ 3-^- - J (48) 



in Landau gauge. To get the contribution from diagram (e), replace the group factor, 3, 
with tan 2 &w 

In this gauge, the numerical value of the contribution to the wave function renormalization 
from diagram (c) is — 0.017 /Pl, while the contribution from the "tadpole" diagram (b) is 
+0.253/ Pl- The tadpole contribution is much bigger, even though diagram (6) vanishes 
in the continuum theory and would not naively be expected to contribute to the wave 
function renormalization at all. The contribution of this diagram almost exactly equals 
the expectation of the "tadpole improvement" argument. 

For a general gauge parameter, where the gauge propagator becomes 

6a + (a - 1)^- 
3 ~ k2 k2 , (49) 

the contribution from (b) changes by a(E/47r)(— 3 +p 2 /2), and the contribution from (c) 
changes by — a(S/47r)(— 3 +p 2 /2) — 3ap 2 £/47r. The mass squared correction is not gauge 
parameter dependent, but the wave function correction is. This will be matched by a similar 
gauge parameter dependence in Zop, so the combination Z^Zq\> is gauge invariant. 

New complications arise in computing the gauge particle self-energy. The value of the 
self-energy must be transverse and rotationally invariant, and it must approach zero at 
p — > 0; but individual diagrams need not separately satisfy these requirements. It is a 
nontrivial check on the calculation if they do, and on the lattice that check will involve the 
use of identities, which are intimately related to gauge invariance. To illustrate some of 
these issues we will present the calculation of the self-energy corrections from scalar loops. 
Two diagrams contribute, diagrams (/) and (g) in Figure |4]. The contribution of (/), which 
vanishes in the MS renormalized continuum theory, is 

7T A:L from (/) = -^-5 AB 5ij f , (50) 

PL J[-k,tt] k 2 



12 



(f) (h) 




(1) 



Figure 4: Contributions to the SU(2) gauge particle self-energy tt a . The curly line is the 
gauge particle, the solid line is the scalar, and the dotted line is the ghost. 



where the 2 arises from the trace over the fundamental representation (in the standard 
continuum notation it would be 1/2, but here the coupling involves the Pauli matrix r rather 
than r/2), and the S AB is in group space. Using cos/cj = 1 — kf/2 and cubic invariance, we 
perform the integrals and get 

2 / E 1\ 

n AyL from (/) = -— <W*j ^— - - J . (51) 

This does not vanish as p — > 0, and in fact does not depend on p at all. It is also independent 
of the gauge parameter a, because no gauge propagators appear in the loop. 
Diagram (g) contributes 

ir A L from (g) = ^8 AB - [ J2fc) i (2fc)^ 

in the lattice regulation and 

^ from {9) = J l 5ab \ L {k-rimk + rW - UV divergence (53) 

in the continuum. The "— UV divergence" subtracts the value at p — 0. Hence the 0(p°) 
contribution is just the lattice value, which is (using (2k)i = 2/cj cos(/cj/2)) 

S AB f 4k i k j cos{k i /2)cos{k j /2) = 5 AB 5 i:j r 4fcf - k\ 
Pl 'l-^f ' " (//-') 2 " Pl A-*,*] 3 (k 2 ) 2 
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The integral involving k\ can be performed using cubic invariance and gives (4/3)(S/47r), but 
it is not immediately obvious how to perform the other integral, or that its value will correctly 
cancel the other 0(p°) parts. The integral is solved using the identity, Eq. (|113|) , which was 
derived from the invariance of J^_ n ^3 ln(k 2 ) on shifting k\ (which can be compensated for 
by a shift in the integration range). It does in fact cancel the other 0(a°) parts arising from 
scalar loops. 

To understand why, note that all 1 loop contributions to the gauge action from scalar 
loops can be computed by performing the integral over the Higgs fields in the Gaussian 
approximation (ie neglecting Al); the scalar loop contribution is ( — l/2)Tr ln(— D 2 + m 2 ). 
This can then be expanded about Tr In —d 2 +m 2 = J,_ n ^3 ln(A; 2 +m 2 ). The required identity 
came from shifting the integration variable in exactly this integral (at m 2 = 0), which is of 
course the same as applying a spatially nonvarying gauge field. In fact, Eq. ( |112j ) contains 
precisely the integrals in Eqs. ( |50|) and ([54]) which contribute at 0(p°). 

Expanding the contribution from diagram (g) to 0(p 2 ) one finds, from Eq. 



5ab 

(3l 



Akikj cos(fcj/2) cos(/cj/2) 



2^2 



v i (h 2 

Y 2fc 2 V2 ' 



1 + 



kik m cos(k/2) cos(fc m /2) 



Irn 



2\2 



The corresponding continuum contribution is 



5 



AB 



Pl k A 



p 2 (p ■ ky 

'2k 2+ k 4 



(55) 



(56) 



which removes the infrared divergence of the lattice version. The difference between the 
lattice and continuum expressions must be of form 



5 



AB 



0L 



{AS^p 2 + Bp^j + C8 i:j p, 



(57) 



simply from cubic invariance. We expect that the answer will be rotationally invariant and 
transverse, B = —A and C — 0, because all gauge invariant, cubic invariant, dimension 4 
operators are; and this will constitute a check on the calculation. The p\ contribution when 
% = j — 1 is 



A+B+C =4 



k 2 {\ - k 2 /A) (k\ - 2 k\{\ - k 2 /A)\ v 

+ , r „,„ 1/ - IR divergence , (58) 



2^2 



Ak 2 



(k 



2\2 



which can be evaluated using Eq. (|115| ) and Eq. ( |120|) in Appendix A; it vanishes. One can 
find A from the p\ contribution when i — j — 1; it equals 



A 



k 2 {l-k 2 /A) (k 2 -2 , k 2 {l-k 2 /AY 



(k 2 ) 2 V 4P (k 2 ) 2 J 

which can be evaluated using Eq. (|116|) and Eq. ( |121|) , and equals 

is u 



IR divergence 



.4 



12 4tt 3 4vr 



(59) 



(60) 
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To evaluate B one takes the p±p2 part when i = 1 and j = 2, which is 

kikl{l-~ki/A){l-~kl/A) 



B = 8 



— IR divergence = —A . 



(61) 



-,-] 3 (P) 4 

evaluating it with the same equations. Hence C = and the result is indeed rotationally 
invariant and transverse. 

There are no new techniques involved in the remaining integrals, so we will simply tab- 
ulate all of the results here. 

Five diagrams, (a) through (e) of Figure |], contribute to the Higgs self-energy. They are 
evaluated above, and equal 

j] 

Pl(k<s>,l - 7r$,c) = _ 6A— from (a) + (62) 



47T 

-S-Sa^+^ + f)^ from(b) + 



E o / . 1 a\ E , , £ 

3a S + " ((-2-2'S + < 9 - 3Q ')S 

/ „ \ S o,l a. E 

6 6 ; 4vr V Un 



(63) 

from (c) + (64) 
tan 2 Q w from (cZ) + (65) 
tan 2 0vk from (e) . (66) 



Seven diagrams, (/) through (Z) of Figure |], contribute to the SU(2) field self-energy. In 
Feynman (a = 1) gauge, they contribute (Sab understood) 



Si 



Sy ( 2 



-2— + - 

4tt 3 



47T 



from (/) + 
+ ($ijP 2 ~ PiPj) 



'IE 1 C 
124vr ~ 34vr 



from (h) + 



4 j 



from (z) 



+ 



(67) 

from ((/) (6*) 

(69) 
(70) 



6 4vr 3 4tt 



6 4^ + 347 



from (j) + 



5y 



-16 



PiPj 



Si, 20 



E 28 

+ y 

8 E 1 
~347 ~ 3 
E 4 



+ Sijpi 



14 E \ 



+ 



1 



47T 



-16 



1 



47T 



17 E ' 



Sid ( 2 



4tt 3 
29 E 
1T4> + 3 
E 1 
'4?r ~ 3 



(71) 



from (jfe) +(72) 



+ 



from (Z) . (73) 
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(m) (n) 

(o) 

Figure 5: Contributions to the U(l) particle self-energy. The wiggly lines are the U(l) 
gauge particle, and the straight lines are the scalars. Diagram (o) only exists in the lattice 
theory. 



Only the last two are gauge parameter dependent. The total is rotationally invariant, trans- 
verse, and free of 0(p°) contributions. 

Only three diagrams, (m) through (o) in Figure [5], contribute to ttb- Scalar loops con- 
tribute the same value as for the SU(2) self-energy, because although the SU(2) contribution 
is enhanced by a group factor of 2, the SU(2) doublet contains 2 fields which are charged 
under U(l). The result is 

'' ^b,l-kb,c) = + from(m) + (74) 



tan 2 Q w 



[l2^- 3 47 J fr ° m(n) + (75) 
■^{SijP 2 - PiPj) from (o) (76) 

The factor of tan 2 Qw arises as follows. Each B field propagator contributes a tan 2 Qw, and 
there is 1 more propagator when a self-energy insertion occurs than when none occurs; in the 
case of diagram (o), there are two extra propagators, but the vertex contributes tan -2 Qw- 
Although this diagram does not exist in the continuum theory, it completely dominates the 
contribution to the self-energy. It arises from the compact nature of the gauge action, and 
its value is given exactly by the tadpole improvement technique at one loop. 

Some 12 diagrams contribute to the scalar four point vertex, in a general gauge. All but 
4 of these vanish in Landau gauge; also the diagrams involving gauge lines can be grouped 
by topology, with the type of gauge particle (A or B) only determining counting factors. 
The diagrams are listed in Figure ^; their contributions are 

Pl(V x , l -V x , c ) = -12A 2 J^ from(p) + (77) 
-(4 + 2(1 + tan 2 9 w f)f + 

Alt 

-(2 + (l + tan 2 w ) 2 )a 2 — from (q) + (78) 

An 

2(2 + (l + tan 2 VK ) 2 )a 2 — from (r) + (79) 

An 
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(t) (t) 

Figure 6: Contributions to the scalar 4-point function, which correct \ L . All but (p) and 
(q) vanish in Landau gauge. The curly lines are SU(2) gauge propagators, the wiggly lines 
are U(l) gauge propagators. 
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-(2 + (l + tan 2 9 w ) 2 )a 2 ^- from (s) + (80) 

(6 + 2tan 2 e,4/)a^A from (t) . (81) 

The a 2 terms all cancel, and the a term absorbs the a term in Z<$, in the computation of 5Xl- 
The Landau gauge value could have been derived by inserting m 2 in the propagators of the 
"tadpole" contributions to the scalar self-energy and finding the 0(m 2 ) term, or equivalently 
by expanding the (tadpole) one loop contribution to the Landau gauge effective potential to 
a sufficient power in m, as discussed in Appendix A. 

Seven diagrams contribute to the scalar- A field 3 point vertex. Again, they can be 
grouped by topology; they are presented in Figure 0. Their contributions are 

Pl(Va,l ~ V A , C ) = (J + tan 2 W ) (-± - |) |j- from (u) + (82) 

1 + tan 2 Q W ) + ^ + from (v) + (83) 



(-l +t an 2 e w )^£-|^j fromH + (84) 

( 3a- — | from (x) . (85) 

V 4tt 6 47r/ v ; v ; 

When using this result to get one must remember that the A field self-energy was only 
evaluated at a = 1; note that there is a nonzero a dependence in Va + tt<i>, and hence also in 

71 A . 

The diagrams contributing to Vb are the same as (u) through (w), but with A and B 
lines switched. The contributions are 

ft,(V B ,L - Vb,o) = (3 + tan 2 e B .)(-5-|)^ from («) + (86) 

(3 + ta„ 2 e w )(-3i- + i| + ||) fromM + (87) 

(3 + tan 2 Q w ) (a— - ——] from (it;) . (88) 

V ' \ ATT 6 47T J v ' v ' 

The total is precisely minus the total 0(p 2 ) contribution to the Higgs self-energy; hence Vb 
and Z$ cancel in the evaluation of Zb and only ttb contributes. The same thing happens in 
the renormalization of continuum, 4 dimensional U(l) gauge theory. 

Finally, there are three contributions to the multiplicative renormalization of the $ 2 
operator insertion, listed in Figure The contributions to L$ are 

(3 L (L^ L -L^ C ) = -6X-f from(y) + (89) 

47T 



a (3 + tan 2 O w ) ^- from (2) , (90) 



which will cancel the a dependence in Z$ when one forms the gauge invariant combination 
Z^Zqp- 

This concludes the evaluation of the relevant diagrams. 
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vanishes vanishes 

Figure 7: Contributions to the SU(2) Higgs vertex. The wiggly lines are U(l) propagators, 
the curly ones are SU(2) propagators, the straight ones are scalar propagators. The last two 
diagrams vanish on integration over the loop momenta. 
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♦ 




♦ 



(y) 



(z) 



(z) 



Figure 8: Contributions to the $ 2 operator renormalization. The dots are $ 2 operators, 
the plain lines are scalar propagators, the curly line is the SU(2) propagator, and the wiggly 
line is the U(l) propagator. The diagrams with gauge particles automatically vanish (at zero 
momentum) in Landau gauge. 

4 Conclusion 

We have argued that the substantial linear in a corrections arising in lattice Monte-Carlo 
investigations of the strength of the electroweak phase transition arise from the difference 
between ultraviolet screening of couplings and wave functions in the continuum and lattice 
theories, and can be cured either with an 0(a) shift in the relation between the continuum 
and lattice couplings and wave functions or equivalently with an 0(a) shift in the interpre- 
tation of lattice data (the values of A and (3l used and the jumps in order parameters). 

To test the technique, we re-examine existing data, the infinite volume extrapolations 
of the jump in the order parameter A(/3l& 2 ) (the difference of $ 2 between phases at the 
equilibrium temperature, extrapolated to the infinite volume limit) for m* H = 35 and 60GeV, 
computed in ||. We applied 0(a) corrections to the data using Eqs. (|2l f )- (|3~0"|) p|. We plot 
the uncorrected and corrected jumps in the order parameter in Figure |9], which also shows 
the two loop perturbative value as a function of A. The corrected data tell a consistent story. 
At m* H = 35GeV the two loop perturbation theory is working fairly well, underestimating 
somewhat the jump in the order parameter, presumably because of transition strengthening 
three loop effects || (though the coarsest lattice seems inconsistent with the other two, but 
see below), and at 60GeV the jump is beginning to fall below the perturbative prediction. 
Before the improvement the points demand a substantial extrapolation; the result of this 
extrapolation is also shown in the figure. Such an extrapolation should successfully remove 
0(a) errors, but it is less optimal than improving the lattice data analytically because the 
extrapolation is numerically expensive due to the demands of the finest lattices, and the 
errors are dominated by the 0(a 2 ) errors in the coarsest lattice (which may also be amplified 
by the extrapolation process) and the statistical errors in the finest lattice. Using a single 
intermediate coarseness lattice and the 0(a) improvements could give smaller statistical and 
systematic errors at less numerical cost. 

For the m* H = 35GeV data the linear extrapolation of the uncorrected data is quite poor, 
with x 2 — 2.9 for 1 degree of freedom. Comparing the corrected data to the perturbative 
estimate, we see that the two finer lattices are following the expected trend of decreasing <fi 
with increasing A; but the data on the coarsest lattice buck the trend, suggesting that the 

3 There is an 0(a 2 ) ambiguity involved in solving these equations-should one compute Za etc. using the 
uncorrected or corrected couplings? We used the corrected couplings, which means that the expressions had 
to be solved iteratively. 
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0.064 0.066 0.068 0.07 0.072 0.074 0.22 0.24 0.26 

lambda (lattice units) lambda (lattice units) 

Figure 9: Jump in /3|$ 2 as a function of scalar self-coupling, before 0(a) corrections 
(squares) and after 0(a) corrections (triangles). The data is that used in ||. The figure at 
the left is the m* H = 35GeV data, at (3l = 8, 12, 20 (from top to bottom before improvement 
and from small A to large A after), and the figure at the right is the m* H = 60GeV data, 
at Pl = 5, 8, 12, 20. In each case the last unimproved datapoint is the result of a linear 
extrapolation of the others. 
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blame for the poor fit falls on the data from the coarsest lattice (/3l = 8). This is probably 
due to 0(a 2 ) effects. At such a small value of A, the phase transition is very strong, and 
the particle masses in the broken phase are quite large. The natural scale of the physics 
involved in the phase transition is quite short, short enough that one should worry about 
corrections from high dimension operators. We can estimate the size of such corrections 
by seeing how much one such effect, a </> 5 //^2 term induced in the effective potential at 
one loop by nonrenormalizable derivative corrections (computed in Appendix A), changes 
the perturbative value. Its contribution to the effective potential, summing only over the 
transverse gauge fields, is 



Adding this to the 2 loop effective potential shifts the minimum (in the lattice units used in 
the figure) from /?|$ 2 = 35.4 to /3|$ 2 = 30.2, which is of the right sign and the right general 
size to account for the discrepancy between the data and the perturbative value. (This is 
not the only 0(a 2 ) effect, so we should not have expected good agreement after correction-it 
can only be used to estimate what accuracy we can expect.) The correction for the next 
datapoint is already smaller by a factor of 0.4. Because of this error, the extrapolation of the 
unimproved data suggest a jump in the order parameter which is larger than the perturbative 
value, by an amount which is inconsistent with the corrected data from the two finer lattices. 
This is an example of the danger of linear extrapolations; 0(a 2 ) systematic errors from the 
coarsest lattice data appear in the final answer, which depends very strongly on the coarse 
lattice data because its statistical errorbars are small. 

The 0(a 2 ) correction from nonrenormalizable operators is smaller for the 60GeV data, 
because a 5 term changes the strength of the phase transition far more for small A than for 
large A; from the one loop perturbative effective potential we would estimate the fractional 
change in the order parameter due to such a term to be oc 0/A ~ 1/A 2 . A small value of A 
demands a very fine lattice, which makes intuitive sense because a stronger phase transition 
involves physics on shorter length scales. Hence, the extrapolation of the 60GeV data much 
better represents the trend in the improved data. 

The technique we have presented here is not restricted to correcting the jump in the 
order parameter $ 2 . One can also correct computed surface tensions, by replacing the naive 
values of Pl and A with the corrected versions in those calculations. Improving the surface 
tension calculation is particularly easy because it only depends on a physical length scale and 
ratios of probabilities for different values of the order parameter, and not on any operator 
insertion which might renormalize. To remove 0(a) errors from the calculation of other order 
parameters, such as the string bit or the Wilson loop, one must still compute the one loop 
corrections to the appropriate operator insertions. 

While we estimate the size of an 0(a 2 ) error above, so far we have said nothing about 
removing 0(a 2 ) errors. These arise not only from the finite renormalization of the bare 
couplings and wave functions, but from tree level nonrenormalizable operators, which modify 
the infrared behavior of the theory at 0(a 2 ) (for instance, the <p 5 term mentioned above). To 
remove them one should start out with an "improved" action rather than the minimal 
Wilson action we discuss here. Unfortunately the ultraviolet effects computed here will be 
different in such an improved action and must be recomputed. 




lattice units , 



or 



6 g 5 a 2 T 
160vr 32 



5 , physical units . 



(91) 
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This still leaves out the 0(a 2 ) corrections to couplings and wave functions which arise 
from two loop graphs. Directly computing these would involve 205 two loop diagrams for 
SU(2) Higgs theory, [] and even more if the U(l) factor is included. It is clearly impractical 
to attempt this renormalization. However, one can absorb the great majority of the 0(a 2 ) 
corrections to wave functions and couplings in a very economical fashion, by using the tadpole 
improvement scheme. First, the action is tadpole improved by the technique developed by 



Lepage and Mackenzie [12] , thereby correcting most of the loop contributions to the wave 
functions and couplings, at every order. Then 0(a) corrections to the wave functions and 
couplings are applied, to compensate for the difference between the full 0(a) corrections and 
those 0(a) corrections which the tadpole scheme will have already made. This removes all 
0(a) errors and most 0(a 2 ) and higher order errors. 

Note finally that while removing all 0(a 2 ) errors may be impossible in practice, removing 
all 0(a 3 ) errors is probably impossible even in principle, because there are 0(a 2 ) corrections 
to the ultraviolet propagators, brought about by their couplings to infrared nonperturbative 
physics, which are not computable and which propagate via the one loop UV corrections 
discussed here into 0(a 3 ) errors in the infrared wave functions and propagators. However, 
from a pragmatic point of view, since the dimensional reduction program itself has a limited 
accuracy, and since it is quite easy in practice to drive 0(a 3 ) corrections below the 1% level, 
this is not of practical importance. 
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A Integrals and Identities 

In this appendix we expand the tadpole integral in powers of mass, evaluating two of the 
resulting integrals numerically. Then we derive a number of identities which relate all other 
integrals needed in this paper to these two. 

The basic integral encountered in finding the effective potential to one loop in the lattice 
regulation with the minimal (Wilson) action is || 

al(rn) = f f * * , (92) 
J[-tt,tt] 3 (In) 6 k 2 + (am) 2 

where, as in the text, k\ = 2sin(fcj/2) and k = Ya=12 3 4 sin 2 (A;j/2). Also, the d 3 k/(2n) 3 will 
be assumed in all integrals in the remainder of the section for notational simplicity. The 
integral has been expanded in a power series about am = through 0((am) 3 ) in 0. We 
repeat that expansion here, but show exactly what integral is responsible for the third and 
fourth coefficients, rather than getting them numerically by evaluating the original integral. 
This is useful because the integral for one coefficient arises in the one loop renormalization 

in Landau gauge, where most of the scalar four point vertex corrections vanish automatically, there 
are 17 diagrams contributing to this vertex, 31 diagrams contributing to the scalar self-energy, 35 diagrams 
contributing to the <I> 2 operator correction, 50 diagrams contributing to the gauge self-energy, and 72 non- 
vanishing diagrams contributing to the gauge-scalar vertex. In a general gauge there are even more. 
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performed in the text, and the other illustrates the influence of nonrenormalizable operators 
on the infrared physics. 

The integral is well behaved about k = and so to lowest order in am <C 1 it is 



1 S 

-tt.tt] 3 k 2 An 



3.175911536 



(93) 



The notation is borrowed from Farakos et. al. || , who found an analytic expression (although 
the value presented above is the result of an accurate numerical integration). This gives the 
(divergent) mass squared correction. 

Adding and subtracting this integral, we get 



al(m) = [ 

All J \ — 7T,7rl' 



(am) 



k 2 (k 2 + (am) 



(94) 



The second integral is infrared dominated, and it is best to add and subtract the analogous 
continuum integral, which can be performed; 



am 



[-n,nf k 2 (k 2 + (am) 2 ) 



[am) 



sr3 k 2 (k 2 + (am) 2 ] 



[am) 



+ 



(95) 



am 



■*,*] 3 \k 2 (k 2 + (am) 2 ) k 2 (k 2 + (amf 

|2 



(am) 



(am) 



S3 k 2 (k 2 + (am) 



» 3 -[-7r,7r] 3 k 2 (k 2 + (am) 2 ) 
am 
An 



(96) 



This gives the famous negative cubic term in the effective potential, and leaves a remainder 
which vanishes as (am) 2 . 

Because k 2 — k 2 = 0(k i ) at small k, the remaining terms are well behaved about k = 0, 
and we can again extract the dominant behavior by adding and subtracting the values at 
am = 0: 



[am) 



am 



am 



k 2 (k 2 + (am) 2 ) k 2 (k 2 + (am) 2 

' ---) 

k* 



[am) 



.(k 2 ) 2 
1 



K 3 -[-7r,7r] 3 k 4 
1 



j?3-[-7r,7r] 3 k 2 (k 2 + (am) 2 ) 
am) 4 x 

1 



-tt^] 3 V (k 2 ) 2 (k 2 + (am) 2 ) k\k 2 + (am) 2 ) J Jw-[-^f k\k 2 + (am) 2 ) 
Following the notation of pf we denote the first integral as 



(97) 
(98) 



[ — 7T,7r]' 
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2\1 



k 4 



SR3_[_ 7rj7r ]3 £.4 



An 



(99) 
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To evaluate it numerically we reduce the second integral to a single integral over a polar 
angle, 

/ JL = _L + J_ cos(#) arctan(cos(0)W0 = 0.0134035706 (100) 

JS 3 -[-7r,7f] 3 K 47T d 7T 4 JO 

and perform the first numerically, directly. The result is £ = 0.152859325. This enters the 
coefficient of an 0(a) contribution to the effective potential which behaves as 4 ; summing 
over species, one can derive V\ in Landau gauge. 

Finally, Eq. (®8p is infrared dominated and 0((am) 3 ), because k 2 — k 2 is order k A at 
small k. To extract this infrared 0((am) 3 ) piece one should expand k 2 to 0(k 4 ), 

£4sin(A,/2) =E*f-^E*i. ( 101 ) 



from which it follows that the contribution — Eq. ( |9~8"D to the effective potential is approxi- 
mately 

2A; 2 + (am) 2 2 \ 




2^2 



32, (102 > 

This generates an 0(a 2 ) term in the effective potential of form 5 and is caused by the (tree 
level) nonrenormalizable derivative terms in the action arising from the choice of lattice 
regulation. Had the expansion of k 2 been free of the 0(k 4 ) part, no such term would appear 
and the first one loop effect beyond the 0(a) correction to A would be an 0(a 3 ) induced 
nonrenormalizable 6 term. (Of course, there would still be corrections arising at two loops, 
for instance an 0(a 2 ) correction to A.) 

The result for the tadpole integral is then 

. . E am £{am) 2 (am) 3 n4n 

and the contribution to the effective potential, / ml(m)dm, is 

S m 2 m 3 fm 4 a a 2 m 5 ^, o R . , 

+ — — + 0(a 3 m 6 ) . (104) 



4vra 2 12vr 16tt 160tt 

The first term is the linearly divergent mass squared correction. The second is the negative 
cubic term which makes the phase transition first order, the third is the contribution to the 
O(a) correction V\, and the fourth is a change induced at one loop from a nonrenormalizable 
derivative correction operator which appears with 0(a 2 ) coefficient at tree level in the lattice 
theory. Of these, the terms with E and £ in them arise from integrals which we need in the 
body of the text. 

Several other integrals arise in the calculations in the text, but a number of identities 
relate them all to the two (numerical) integrals, Eq. (|9~3D and Eq. (|99|). Some of these 
identities involve trigonometric manipulations of numerators, for instance using cos 2 (fcj/2) = 
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1 — kf/4 and cos(A;j) = 1 — kf/2; we will not bother to list these here. Some use cubic 
invariance; 



L 



k{ ir kf + ki + ki if 1 is 



n,nf (k 2 ) 2 3 7[-^] 3 (k 2 ) 2 3J[-^fk 2 3 4vr 



Similarly, 



and 



f S = i 

J[-^f k 2 3 



(105) 



(106) 



L 



k 2 1 £ 

1 continuum version = -— . (107) 



-7r,7r] 3 (A; 2 ) 3 3 4vr 

Here "- continuum version" means that kf/k 6 should be subtracted from the integrand and 
the integral of kf/k 6 outside the first Brillion zone should also be subtracted. In what follows 
the meaning will be analogous. 

Cubic invariance cannot be used to evaluate the integral over kf/(k 2 ) 3 . To do so we 
must take advantage of the periodicity of the integrand on In shifts in any fcj. For e an 
infinitesimal, the periodicity of k\ ensures that 



/- 



1 



[-n,nf k 2 Jh^f (2 sin 2 ((A; 1 + e)/2)) 2 + k\ + k\ 



(108) 



Now, excising a small neighborhood about the origin to avoid the singularity there, one can 
expand in e. At order e 2 , one finds after some work that 

= / — = — + — ^— + 4^ J— + surface term , 109 

J[-n,n] 3 -excision ^ (k 2 ) 2 2 (jfe2)3 (jfe2)3 (£2)3 J , \ J 

where the surface term arises on the boundary of the excised region. It cancels the linear 
infrared divergence occurring from the two terms with 0(1/ k A ) infrared behavior, so what 
remains is the difference between these terms and their continuum analogs, eg Jl/(k 2 ) 2 
becomes Eq. (|99p . Using the previous integrals and identities we can integrate all but the 
last term, so we find that 

r kf 1 E 1 £ . . 

/ = + -— . 110 

J[-n,nf (k 2 ) 3 6 4vr 3 4tt v 1 

Because k\ + k 2 + k 2 = k , it follows from this integral, Eq. ( |105 ), and cubic invariance that 

kfk 2 IS 1 £ 



-tt.w] 3 (jfe a )3 12 4vr 6 4vr 



111) 



Repeating the above technique, but starting with the integral over ln(k ) rather than 
l//c 2 , one finds that 

r -2k 2 + k 2 k 2 + Ak 2 - fcf t s 

0= / J 1 112 

Jl-nri 3 (k 2 ) 2 



26 



or 

f iL = i_?s, (113 ) 

J[-7r,n] s (k 2 ) 2 3 3 4vr' v ' 

and hence also 

A--,-] 3 (k 2 ) 2 3 4tt ( j 

Continuing the expansion to fourth order in e, and using the previously derived identities, 
one obtains after considerable algebra that 

kUl - k 2 /A) 2 

■ — continuum version— 



(P) 4 

2,1 l2/n\n L2 



f jfc?(l-jfc?/2)(l-fc?/4) . . 1 IS 1 C ,„rs 
/ = continuum version = 1 . 115 

Jl-n^f (A; 2 ) 3 48 8 4tt 8 4tt v 1 

Shifting the k 2 integral by 5, the condition that the 0(e 2 5 2 ) term should vanish gives that 
kfk 2 (l-k 2 /4)(l-~k 2 /4) 

= continuum version— 

b^] 3 {k 2 Y 

r k 2 (l - fc|/2)(l - k 2 /A) IS 1 i . . 

I ~ continuum version = — . (116) 

J[-K,n] 3 (k 2 ) 3 32 4tt 8 4vr v ; 

By shifting ki by e in the integral 



^1 

l[-n ,^] 3 k 2 

and expanding to second order in e, one finds using Eq. ( |114j ) that 

[-^] a (jfc2)3 ~ 6 47 ~ 3 47 ' 
and hence 

k\ 1 S 4 £ 



117) 



'118) 



(119) 



/[-^] 3 (k 2 ) 3 3 4tt 3 4tt 

by applying the same technique used to get Eq. ( |111| ) and using Eq. ( |113| ). With these 
integrals we can derive that 

r kj{l-k 2 /2){l-k 2 /A) 1 1 S 1 £ , . 

/ = continuum version = 1 — . (120) 

Jl-n^f (k 2 ) 3 24 4 4vr 4 4vr' v ; 



and 



f k 2 (l - k 2 /2)(l - fei/4) . . IS 1 e 

/ tl—L - continuum version = 1 — . (121) 

Jl-n^f (k 2 ) 3 16 4vr 4 4tt v ' 

Hence we know the value of each integral in Eqs. ( |115| ) and ( |116| ) separately. 
No other identities are needed for to complete the 1 loop renormalization. 



27 



B Feynman rules 



Most of the Feynman rules needed in the calculation appear in [TJ|; the measure insertion 



(diagram (h) in Section [3]) is smaller by a factor of 2/3 from the one printed there, which 
is for SU(3) rather than SU(2), and similarly in the four point vertex in Eq. (14.44) one 
must set (Iabc to and the (2/3) in front of Sab^bc + ^ac^bd + • • • should become 1. Our 
conventions for the scale and the gauge wave function mean a should be replaced by 1 and 
g by 2 everywhere. Also, the sign for the ghost coupling to two gauge particles is backwards 
there. 

What remains are the scalar propagator and the vertices involving the scalar. The prop- 
agator, in the ultraviolet where m is negligible, is just 5 a b/p 2 - The 4-point vertex is 



— 2X(5 ab 5 cd + 5 ac 5 bd + 8 ad 8 bc j 



122) 



where a, b, c, d are the species types of the incoming lines. The vertices involving gauge 
particles, and the pure U(l) four point vertex, are given in Figures to E| The notation 
is that {i) a b means the action of i on the species type, eg species type 1 becomes type 2 and 
type 2 becomes — type 1, etc. Any % by itself is the imaginary number which arises when 
one Fourier transforms an odd function. The vertices in Figure ^ do not have continuum 
analogs. Their role in generating extra tadpole diagrams is quite important on the lattice; 
for instance they dominate the corrections Vb to the U(l)- Higgs vertex. 
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P,i qj 




r,k s,l 
r,i,A 



£(p5 - p5 )(q8 - qS ) 

Z m n m in n im m J n n J 1 " 

' (r8 - rS )(s 8 - sS ) 

v m kn n km x m In n to' 



-(i t a ) i p-q 5(p+q+r) 

ab i i 



p,a 



q,b 



r,i 



-(i) ab i P ; - q ; 8(p+q+r) 



i i 



p,a 



q,b 



Figure 10: Feynman rules for the U(l) four point function and the three point gauge-Higgs 
couplings. 
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r,i,A 



p,a 



q^b -2 5 AB 8ij 5 ab cos( ^ i ) 5 (p+q+r+s) 



p,a 



q,b 



-2 5 ab 5.. cos( ^lJ 5 (p+q+r+s) 



p,a 



q,b 



2 (i)ac TA ) C b C0S ( ^ *) § (P+q+r+s) 



Figure 11: Feynman rules for the four point gauge-Higgs couplings. 
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r,i,A 



p,a 





s,j,B 



t,k,C 



q ,b -y- ( 5 AB 5 CD + 5 AC 5 BD + 5 AD 5 BC ) 5 ijk 
(ix D ) ab i P : - q 5 5 (p+q+r+s+t) 



i i 




^ Sab 5 ijk (i) ab i p. - q . 5 (p+q+r+s+t) 



s,j,A 



t,k,B 



r,i,A 



p,a 



^ (i TA )ab 6 ijk 1 R -q ; S (p+q+r+s+t) 



so 



t,k 



r,i 



i i 



p,a 



q,b (i) h 5 iik ip ; -q ; 5 (p+q+r+s+t) 



'ab ijk *i 

s,j t,k 

Figure 12: Feynman rules for the five point gauge-Higgs couplings. 
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C Results in other theories 



Here we present results for the renormalization of the abelian Higgs model with N scalars, and 
for SU(2)xU(l) theory with adjoint scalars (the Weinberg-Salaam model after dimensional 
reduction but before integrating out heavy modes). 

The action of the 3-D abelian Higgs model on a lattice is 



E 



N 



Z B J2(l-Re Pij ) + Z^H 

i<j 6=1 i 



- $ ab (x)u iac (x)§ cb (x + + 



rn 



HL , 



EC 

6 



+ zl 



t(E^) 2 



(123) 



Here a and c are indicies over the real and imaginary part, which are implicitly summed 
over when $ 2 appears, and b is an index over the N scalar species, which are assumed to 
have the same mass squared and an SU(N) symmetric quartic interaction. The finite a 
renormalization only involves diagrams computed already in the text; the results are 



Pl(Z b - 1) 
f3L$m 2 HL 



Pl(Z< 



+ / 6 A L -(8 + 2iV)Ai-2 
6 An v 

3 \24 4tt QAtvJ ' 
(2 + (2 + 2iV)A L ))^, 



An 



(124) 
(125) 

(126) 
(127) 



Next we consider SU(2)xU(l) theory, with an adjoint SU(2) scalar Aq and an adjoint 
U(l) scalar Bq. If one is interested in very high precision calculations of the electroweak phase 
transition strength one must include these, because the integration over these fields (the so 
called heavy modes) is lower precision than the integration over the nonzero Matsubara 
frequencies (the superheavy modes). 

The integration over the nonzero Matsubara frequencies generates mass terms for these 
modes. Interaction terms already exist at tree level, and new ones are generated at one loop 
in the dimensional reduction. The new mass terms will be denoted m\ and m 2 B , and the 
new interaction terms in the Lagrangian are 



^(Alf + >f(BlY + h -feAl + tf&Bl + h -fAlBl 



-A) • r) Qb $ 6 



'128) 



In the minimal standard model, at lowest order in aw, the bare values for the new couplings, 
in lattice units, are [|J 

(129) 
(130) 
(131) 
(132) 
(133) 
(134) 



Aal 


= o, 


Abl 


= o, 


flAL 


= 2, 


h B L 


= 2 tan 2 ®w 


ABL 


= o, 


K 


= 2 tan Qw , 
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but in what follows we will treat the general problem in which their values are arbitrary. 
Computing the 0(a) renormalizations involves no topologically new diagrams, only combi- 
natorics. The couplings and wave functions presented in Section |2| are modified by 



3,2 



1 s 


4 t 




3 4vr 


\h% L - 




3 L 




- 2 h AL + 


\h B L 







added to (3 L (Z A -1), (135) 

, ,w. i u, ., , added to (3 L 5X L , (136) 

1 ddedto (3 L Sm 2 HL . (137) 

The corrections to Z B and Z$ vanish. 

We denote the new wave function corrections as Z A0 and Z B0 ; the notations for the 
coupling corrections are obvious. The new couplings and wave functions renormalize by 

P L {Z M -1) = 24^ + ^-5A^, (138) 

/3 L (Z B0 -1) = "3A BL |^, (139) 

P L Sm AL = (l6 + 5\ L + 2h AL + ~h AB A—, (140) 

/3 L 5m 2 BL = ^3\ BL + 2h BL + h ABI )j^-, (141) 

PlSXal = - (32 + 11X AL + h\ L + h 2 ABI )j |- + X AL (48^ + ~J , (142) 

PlSXbl = -{lOX 2 BL + hl L + h 2 ABI )jj-, (143) 
P L Sh AL = -(8 + 6X L h AL + 5X AL h AL + 2h 2 AL + hi BL h ABL + 2ti£) J- + 

h AL ^(33 + 3 tan 2 8^)^ + ( y + ~ tan 2 6 w )^j , (144) 
h5h BL = - (6X L h BL + 3X BL h BL + 2h 2 BL + ^h AL h ABL + 2h'lj + 

h BL f( 9 + 3tan 2 e w ,)^ + ( ^) — j , (145) 

^L^h ABL = - (5X AL h A BL + 3X BL h ABL + 2h ABL + 2h AL h BL + 2hf) + 

(3 L 5ti L = h' L ( (-2X L - h AL - h BL - h ABL + 21 + 3 tan 2 Q w ) -L+ 



47T 



7 + tan 2 Q w £ ' 
6 4tt 



(147) 



Here the wave functions already include corrections to make the operator insertions have 
their naive normalization. In fact the operator insertions will be mixed by a matrix with 
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0(a) off diagonal terms; what we use above are the diagonal components of the matrix. 
The off diagonal terms would only be important if the A or B fields took on significant 
condensates at the phase transition, which they should not, because in the realistic case 
they should be given substantial (Debye) masses. They would also be important if we were 
interested in the jumps in these condensates, but this is also not of interest. We have not 
computed the off diagonal terms here. 
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